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AERONAUTIC SYMBOLS 


I. FUNDAMENTAL AND DERIVED UNITS 



Symbol 

Metric -W 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbrevia- 

tion 

Length 

■ l 

t 

F 

meter 

second 

m 

s 

kg \ 

foot (or mile) 

second (or hour) 

ft (or mi) 
sec (or hr) 
lb 

Force 

weight of 1 kilogram 

weight of 1 pound 


P 

V 

horsepower (metric) 


horsepower 

hp 

mph 

fps 

Speed 

/kilometers per hour 

(meters per second 

kph 

mps 

miles per hour 

feet per second. __ _ 


2. GENERAL SYMBOLS 

Weight==m<7 v Kinematic viscosity 

Standard acceleration of gravity =9. 8066 5-m/s 2 p Density (mass per unit volume) 

or 32.1740 ft/sec 2 Standard density of dry air, 0.12497 kg-m~ 4 -s 2 at 15° C 

,, W - and 760 mm; or 0.002378 lb-ff 4 sec 2 

mass— g Specific, weight of “standard” air, 1.2255 kg/m 3 or 

Moment of inertia=m/c 2 . (Indicate axis of 0.07651 lb/cu ft 
radius of gyration k by proper subscript.) 

Coefficient of viscosity . 

3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
•Span 
Chord 

Aspect ratio, g- 
True air speed 
Dynamic pressure, ip I 72 

Lift, absolute coefficient 

Drag, absolute coefficient 

Profile drag, absolute coefficient Dn Q =~^ 

Induced drag, absolute coefficient 0 Df = ^ 

D b 


Parasite drag, absolute coefficient Cdj> = ^§ 


Cross-wind force, absolute coefficient Ca— 


C 

a s 


it 

Q 

o 

R 


a 

e 

a 0 

ott 

Ota 


Angle of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) - ' 

Resultant moment 
Resultant angular velocity 

Reynolds number, p — where l is a linear dimen- 


sion (e.g., for an airfoil of 1.0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) ' . 

Flight-path angle 
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FREQUENCY-RESPONSE METHOD FOR DETERMINATION OF DYNAMIC STABILITY 
CHARACTERISTICS OF AIRPLANES WITH AUTOMATIC CONTROLS 

By Harry Greenberg 


SUMMARY 

A frequency-response method for determining the critical 
control-gearing and hunting oscillations of airplanes with 
automatic pilots is presented. The method is graphical and 
has several advantages over the standard numerical procedure 
based on Routh’s discriminant. The chief advantage of the 
method is that direct use can be made of the measured response 
characteristics of the automatic pilot. This feature is especially 
useful in determining the existence, amplitude, and frequency 
of the hunting oscillations that may be present when the. auto- 
matic pilot has nonlinear dynamic characteristics. 

Several examples are worked out to illustrate the application 
of the frequency-response method in determining the effect of 
automatic-pilot lag or lead on critical control gearing and in 
determining the amplitude and frequency of hunting. It is 
shown that the method may be applied to the case of a control 
geared to airplane motions about two axes. 

INTRODUCTION 

The increased use of automatic control on aircraft (espe- 
cially on pilotless aircraft) lias focused attention on the 
proper design of the aircraft and associated control systems 
with a view toward obtaining satisfactory dynamic stability. 
The control system (autopilot) consists of the gyro, phase- 
shifting device, and servomotor. Factors in the control 
system that determine stability characteristics are control 
gearing, lag in the servomotor, and lead in the phase- 
shifting device. The stability characteristics of the air- 
plane depend on airplane configuration and mass distribution. 

The purpose of the present report is to give a method for 
analyzing the dynamic stability of an airplane with auto- 
matic control which makes direct use of the observed dynamic 
characteristics of the automatic pilot (herein called auto- 
pilot). The method separates the characteristics of the 
autopilot from those of the airplane. It, therefore, easily 
reveals the effects of modifications to the autopilot such as 
adding lead. The procedure is largely graphical. The 
method is similar in certain respects to that of Nyquist 
(reference 1), which was devised for electronic circuits but 
which has also been applied by some workers to the design 
of servomotors. 

The frequency-response method was previously applied by 


Jones (reference 2) to calculate the hunting produced in 
airplanes with “flicker,” that is, on-off control. If the servo- 
motor (herein called servo) has a linear lag characteristic 
(lag independent of amplitude), then there will usually be an 
upper limit to the control gearing above which the airplane 
will be unstable. The frequency-response method deter- 
mines this critical control gearing and the corresponding 
frequency for any type of lag or lead in the autopilot. It 
also indicates the changes to be made in the servo which will 
improve the stability of the airplane. 

The utility of the frequency-response method is most ap- 
parent in the case of an actual servo with nonlinear lag 
characteristics. This method allows the use of measured 
characteristics of the servo that might be inconvenient to 
represent by a mathematical formula. 

TERMINOLOGY AND SYMBOLS 

The word “lead” is used in this report in two ways. It is 
the phase angle of control-surface deflection 5 ahead of air- 
plane deflection, say angle of pitch 6, and it is also used to 
indicate a device that causes an increase in the phase angle 
of lead of S ahead of 8. These phase-shifting devices are 
herein called first-derivative lead and second-derivative lead. 

The expression “linear autopilot” is used to indicate a servo 
that is acted upon by a force proportional to the input signal 
(airplane deflection) and resisted by a force proportional to 
the displacement and velocity of the output (control-surface 
motion) . 

The function of the automatic pilot is to apply a corrective 
control deflection in response to any deflection of the air- 
plane. The “control gearing” is the ratio of the applied con- 
trol deflection to the airplane deflection for very slow 
deflections (static condition). There may exist a “critical 
control gearing,” which in nearly all practical cases is an 
upper limit beyond which dynamic instability occurs. 

The following symbols are used: 
a ratio of control deflection due to second derivative of 
airplane displacement to that due to airplane 
displacement 

r ratio of control deflection due to first derivative of 
airplane displacement to that due to airplane 
displacement 
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R amplitude of airplane oscillation produced by unit 
amplitude oscillation of control surface 
5 r amplitude of control-surface oscillation required to 
maintain unit amplitude of airplane oscillations (1 /R) 
d p amplitude of control-surface oscillation produced by 
autopilot in response to airplane oscillation, divided 
by amplitude of airplane oscillation 
S p b v for unit control gearing 
k„ critical value of control gearing 
k\ spring stiffness factor of servo and control surface 

ko factor relating airplane deflection and force on servo 

k control gearing, ratio of 5 to 6 in static condition 

(k 2 /ki) 

c viscous damping factor of servo and control surface 

l viscous lag factor which depends on c and ki (c/ki) 

m mass of movable part of servo and control surface 

t time, seconds 

T period of oscillation, seconds 

5 angle of control-surface deflection in direction to 

reduce angle of pitch of airplane 
phase angle of lead of 8 ahead of 6 when oscillating 
airplane forces control surface to oscillate 
e r phase angle of lead of 5 ahead of 6 when oscillating 
control surface forces airplane to oscillate 

6 angle of pitch of airplane 

0 signal fed into servo 

cj> angle of bank of airplane 

ip angle of yaw of airplane 

<l> gyro displacement produced by airplane deflection 

o) angular frequency, radians per second (2irXFrequency) 

a n angular frequency of servo and airplane when k=k cr , 
radians per second 

w„ natural angular frequency of servo, radians per second 
r angle of tilt of gyro axis from X-axis of airplane 
D differential operator (d/dt) 

Subscript: 
max maximum 

BASIC PRINCIPLES OF FREQUENCY-RESPONSE METHOD 


gives the amplitude of 5 required to sustain unit amplitude 
oscillation in 6 and is denoted by <5 r ; that is, if 

5=5 t sin at 

then the response will be 

0=sin (c of— t r ) 

A plot of 8 r and e T against a is then made and is referred to 
as the frequency response of the airplane. 

The frequency response of the autopilot is obtained by 
measuring or computing the response of the control surface 
to a steady oscillation of the airplane. The response may 
be calculated if the dynamic constants of the autopilot are 
known, but it is usually easier to measure the response by 
oscillating the autopilot in the laboratory. These calcula- 
tions or measurements give the response to 

0 = sin at 

and this response is expressed by 

5 = sin (cof-f- 

wlieie is the amplitude of 5 for unit control gearing and 
c„ is the phase angle of lead of 5 ahead of 6 when the airplane 
is forcing the control surface to oscillate. 

If at some value of the angular frequency a the condition 

e p =e r 

exists, then the airplane connected to the autopilot will 
oscillate continuously at this frequency at constant ampli- 
tude, provided the control gearing is made equal to S T /S Pi at 
this value of a>; that is, the critical control gearing is given by 



at co where e T =e p . This condition of neutral stability is 
shown schematically in figure 1. An example of a solution 
for k CT is shown in figure 2. In this case co„=6.2 and 


Jr 

Sp 1 


at this value of co; therefore k c 


0.33 
0.91 = 


= 0.36. 


The response of the airplane to a sinusoidal control motion 
may be measured in flight or may be computed from the 
equations of motion of the airplane. (See appendix A.) 
The airplane response is sinusoidal and of the same frequency 
as the control motion and has an amplitude and phase that 
depends on this frequency. If the control motion is given by 

8 = sin at 

then the response in pitch, for example, may be expressed as 
d=R sin (at— e T ) 

where R and e T both depend on a, and R is the ratio of the 
amplitude of 6 to the amplitude of 8 when the control surface 
is forcing the airplane to oscillate. Changing the amplitude 
of 8 produces a proportionate change in the amplitude of 6 
but does not affect the phase angle e r . The reciprocal of R 


Assumed 



Figure 1 . — Response of airplane to sinusoidal control motion and response of control to 
sinusoidal airplane motion. If e p = t r and fc=/c cr =~> the autopilot will sustain the 

d P j 

oscillations. 
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Figure 2.— Illustration of frequency-response method of determining critical control gearing. 
ajer=C\2; ft er— 0.36. (.Dotted line indicates where « p =«r.) 


In almost all cases control gearings greater than k CT result 
in instability and control gearings less than k CT provide 
stability. Thus, k CT is usually an upper limit to the permis- 
sible control gearing. More generally, the stable side of 
the boundary is determined by the relative slopes of the e T 

and curves at their intersection. If ^ is greater than 


the system is stable below k CT and vice versa. A proof 

of this rule is given in appendix B. 

In the preparation of figure 2 and the examples that follow 
(figs. 3 to 8) a particular airplane configuration was assumed 
for the computation of 8 r and e r . (See appendix A for 
methods.) This hypothetical pilotless airplane has four 
fins 90° apart, a horizontal pair for pull-ups and a vertical 
pair for turns. The fins constitute the only lifting surfaces 
of the airplane and are equipped with ailerons connected to 
a roll gyro and servo that maintain the airplane in a position 
in which one pair of fins is always horizontal. The pitch 
and yaw control surfaces are trailing-edge flaps on the fins 
and are connected in pairs, one pair on the horizontal fins 
connected to a pitch gyro and servo and the other pair on 
the vertical fins connected to a yaw gyro and servo. It is 
assumed that the pitch and yaw motions do not interact 
with each other or with any accidental rolling motion that 
826251 — 19 2 



Figure 3.— Effect of lag in the servo on the frequency response of a critically damped linear 

2 

servo. Natural frequency of servo «„=»y 

may exist. Figures 2 to 6 show the elevator-pitch stability. 
Rudder-yaw stability may be expected to be very nearly 
the same as elevator-pitch stability because of the configura- 
tion of the fins. 


LINEAR AUTOPILOT 

The motion of a control surface actuated by a simple linear 
autopilot is similar to that of a mass-spring-dashpot and may 
be represented by the equation 

7ftD 2 8~{-cD5-\-ki5—k26 (1) 

where each of the factors m, c, and k x is the sum of two parts, 
one part due to the servo and one part due to the control 
surface. For example, k x is the sum of the spring constant 
of the servo and the aerodynamic hinge-moment constant 
of the control surface. The term k 2 0 represents a force 
proportional to the deviation in angle of pitch d applied 
without lag. 

Actually, the effect of control-cable flexibility is to intro- 
duce an additional degree of freedom into the system. This 
effect is probably small and may be taken into account 
either analytically or by flight measurements as mentioned 
in the section “Nonlinear Autopilot.” 
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of a critically damped linear servo. Natural frequency of servo a>„*= y 


Equation (1) may be expressed as 

5 0+ F I °+S D ’)=ll» O) 

where k 2 /k i is the control gearing k. The natural frequency 
of the servo is given by 


c 

and the damping depends on the quantity l=jr' With these 
substitutions, equation (2) becomes 

s ( i + io +5)= t » 

The term ID represents viscous damping and the term 
D 2 / represents inertia reaction. 

When the autopilot is oscillated in pitch, the control 
response 5 lags behind 6. The lag usually decreases as the 
lag factor l decreases and as the natural frequency u n 



Figure 5. — Effect of lag and lead on the critical control gearing with a critically damped 


linear servo. Airplane characteristics of figure 2. 


1+r.P 


i+l-n+fi' 


i=- 


increases. It is convenient to choose equal to 2/1, a value 
that gives critical damping. With this choice, the phase 
angle e P and the relative amplitude for unit control gearing 
<hj of the response are presented in figure 3 for four values 
of l. Values of S r and e, for the assumed airplane are also 
given so that k cr can be obtained for each value of l. 

If first-derivative lead is added to the system, the equation 
of motion relating 5 and d becomes 

8 (l+lD+~^)=k(l+rD)d 

Figure 4 shows S Pl and e„ for the same values of l and u„ as 

figure 3, but with r— g- Comparison shows that the first- 

derivative lead produces an angle of lead that approaches 
90° at high values of the angular frequency co. This amount 
of lead is enough to override the lag for moderate values of l 
at low frequencies but not at high frequencies. 

Values of k CT obtained from figures 3 and 4 are shown in 
figure 5 plotted against the viscous lag factor l. Considerable 
increase in the stable range of control gearing is evident for 
small values of l with the first-derivative lead. Note how 
this improvement in stability diminishes at large values of l. 

If both first-derivative lead and second-derivative lead 
are added, the response of the control is given by the equation 

S=k 1+rD +ag o 

1+lD + £ 
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Figure 6. — Effect of first-derivative lead and second-derivative lead on the lag of the servo 
as a function of frequency. 


It can be shown that if 
and 

W> \ (3) 

lead will exist at all frequencies. In this situation the air- 
plane cannot oscillate regardless of the value of k or, stated 
differently, the critical control gearing is infinite. Figure 6 
shows the resultant lag or lead for a critically damped servo 
for which co„=20 radians per second and for various values 
of r and a. The top curve is for values of r and a that satisfy 
expression (3) and, therefore, result in lead over the entire 
frequency range. 

NONLINEAR AUTOPILOT 

Even if the servo is constructed to give a response pro- 
portional to a steady disturbance, its response to an oscillating 
disturbance is not proportional in practice to the amplitude 
of the disturbance because of the nonlinear dynamic character- 
istics of the servo. The lag is then a function of both am- 
plitude and frequency. By finding experimentally the 
frequency response of the autopilot for a number of amplitudes 
and for a given control gearing, it is possible to determine the 
amplitude and frequency at which the airplane will oscillate 
when coupled to the autopilot at that particular value of 
control gearing. The condition for steady oscillations 
(hunting) is that § p = S r and c p = e r at some frequency. The 
frequency and amplitude of the steady oscillations of the 
airplane are the values of co and 6 max at which o p =S r and 
€ P — 6 r . 



Figure 7 - Effect of amplitude on the frequency response of a servo, as measured on an 

oscillating table. 


The stability of these constant-amplitude oscillations has a 
slightly different significance from that of linear systems. 
In linear systems the oscillations are, in general, either 
damped or undamped regardless of the amplitude. In the 
case of the hunting that may exist with a nonlinear servo 
the oscillations are said to be stable if, after a disturbance 
from their steady value, they tend to return to that steady 
value. The criterion for the stability of the steady oscilla- 
tions follows from that for the linear autopilot and is 

d — 

d(e T —e P ) 8 r 
d(*l dO,;, a X 

The solution obtained is not exact because the autopilot 
response to sinusoidal airplane motion is not a pure sine wave 
but is distorted; however, the existence and approximate 
amplitude and frequency of hunting oscillations of the funda- 
mental can be determined by the frequency-response method. 
The values of 8 P and e p used when the response is not sinusoi- 
dal are the values of the equivalent sinusoidal reponse which 
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Figure 8.— Plot of 5 P /d r and against dmax obtained from figure 7, showing the amplitude 

and frequency of the hunting. 

supplies the same energy per cycle and the same impulse 
per half cycle as the observed response. This condition leads 
to the following expressions for t P and o v : 


yA*+B* 

v Omar 

. _! B 

€j,= tan 

where 



The value of S in these integrals is the measured control- 
surface deflection in response to an oscillating motion of the 
airplane represented by 9—9 max sin wt. It turns out that B 


is exactly equal to the cosine coefficient of the fundamental 
frequency component of the Fourier series approximation 
to 8; however, A is not in general exactly equal to the first 
sine coefficient of the Fourier series. Ordinarily the oscilla- 
tion in pitch is performed in the laboratory by mounting the 
autopilot on a table that can be made to oscillate. 

A more exact method of determining the autopilot charac- 
teristics, which includes the effects of control-system flex- 
ibility and hinge moment, is to measure the response of the 
control surface and airplane hi flight to a sinusoidal input 
signal to the servo. In order to obtain the frequency response 
of the control system to sinusoidal input signal, it is necessary 
to divide the vector output 8 by the vector sum of input 
signal due to 9 and input signal applied directly. Thus, the 
dynamic characteristics of the airplane and control system 
can be obtained at the same time. If the aircraft is not 
available for flight tests, one can “synthesize” the complete 
servo-control-system response from ground measurements of 
the servo and estimates of control-system flexibility, inertia, 
and aerodynamic hinge moments. 

A typical application of the method for determining ampli- 
tude and frequency of hunting of an aircraft and autopilot 
with stops in the control system is shown in figures 7 and 8. 
The airplane used in this example is slightly different from 
that considered previously and, therefore, the 8 r and e r curves 
are slightly different from those in figures 2 to 6. Figure 7 
contains the phase and amplitude of the experimental 
response of the autopilot plotted against frequency for a 
series of amplitudes of the input and also the corresponding 
calculated curves for the aircraft. For each amplitude the 
values of 8 P and S r are obtained at the frequency where 
£„=e r . The ratio S P /8 r and the value of co cr (where e P =e r ) are 
plotted against amplitude of pitch 9 max in figure 8. The ampli- 

tude and frequency of hunting occur where ^==1. In this 

example the amplitude of the hunting is about 2.5° and the 
frequency is 18.7 radians per second. 

Nonlinearity in the autopilot usually causes an increase of 
lag with amplitude. This increase may result in a reversal of 
the favorable effect of first-derivative lead as the lead factor 
r is indefinitely increased and suggests an optimum value for 
r, a conclusion that is supported by flight experience. This 
result may be traced to the counteracting effects of lead and 
amplitude introduced by the first-derivative lead. The first- 
derivative lead shifts the phase of the signal fed into the servo 
so that the signal leads 6 by an amount tan 1 no and increases 
the amplitude of the signal by the factor -yJl+rW. The 
change of lead with r may be stated symbolically as follows : 
If the dependence of e„ on 6 and w without first-derivative 
lead is expressed by 

ej ,=j{9, eo) 

then the addition of first-derivative lead changes the lead to 
e, >=/(©, “) + tan _I rco (4) 
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where 

0=0 VF +rV 

Differentiating equation (4) with respect to r gives 

dtp de P d / j— r- -g -A | 

dr dd dr \ ^ 03 ) 1 +r 2 c d 2 


_de p e^r g 

Vl+rW^l+rV w 

If r is very small, the first term on the right-hand side of 
equation (5) is negligible, so that ^ is positive. The addi- 
tion, therefore, of a small amount of first-derivative lead 
results in a change in e„ in the direction of more lead. A 

similar analysis of the effect of r on 8 P shows that < ~^f =0 

when r=0. For small values of r, therefore, the effect of 
adding first-derivative lead is the same as for the linear servo: 
an improvement in dynamic stability and a reduction in the 
amplitude of hunting is to be expected. For large values of r, 
however, the lag due to increase in amplitude (first term) 
tends to counteract the direct lead (second term) and may 
cause more violent hunting. Some optimum value of r may 
therefore exist. A direct determination of its value may be 
made from the autopilot for various values of r, 0, and u>. 


APPLICATION OF METHOD TO CONTROLS GEARED 
TO MOTION ABOUT TWO AXES 

Sometimes the aileron or rudder is geared to both the angle 
of bank and the angle of yaw. Usually this gearing is ac- 
complished by tilting the gyro in such a way that the rotation 
of the gyro is affected by bank and yaw according to the 
formula 

$ = <£ cos r— [— sin T 


In order to apply the frequency-response method in this case 
it is necessary only to calculate the response in bank and 
yaw of the airplane to a sinusoidal control motion and to 
combine them according to the preceding formula to obtain F. 
This formula yields the 8 r and e T curves. The other pair 
of curves required — the 8 P and t p curves — are obtained in 
exactly the same way as before, namely, by oscillating the 
gyro about its sensitive axis. The critical control gearing so 
determined will be the critical ratio between 5 and F. 

CONCLUDING REMARKS 

The frequency-response method of analysis is a useful 
graphical means of determining oscillation characteristics of 
an airplane equipped with an automatic pilot. If the servo 
is linear (ideal case), the critical control gearing beyond 
which increasing oscillations take place can be readily de- 
termined. If the servo has nonlinear characteristics (prac- 
tical case), the existence, amplitude, and frequency of steady 
hunting oscillations can be determined approximately from 
the measured frequency response of the servo and the com- 
puted frequency response of the airplane. 

In general, the use of first-derivative lead (a phase-shifting 
device) has a favorable effect on the dynamic stability. 
Large amounts of first-derivative, lead or phase shift may be 
destabilizing due to the increase in the signal amplitude pro- 
duced by this phase-shifting device. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., December 9, 1946. 


APPENDIX A 


CALCULATION OF RESPONSE OF AN AIRPLANE TO SINUSOIDAL CONTROL MOTION 


The response of the airplane in pitch to sinusoidal elevator 
motion will be used as an example. The basic equations are 
derived in reference 3 but with somewhat different notation 
and without the effects of control deflection, rate of change 
of angle of attack, and power. The equations of motion, 
for a given control deflection, are 

(x u -\- TD)u-\-x a a-\-Xg6=0 

z u u+(z a + TD)a+(z e —TD)d=—z s 5 

c m u + {C\ + C„ 1d D) a + (c mm D - 2Y k* D 2 ) 0= - C ms 5 


where 


_ n 1 (IT 

**-*'•> pVSdV 

_ m 

T ~pVS 

(O h —C D ^j 


x e = 


Or 


Zu = C L 

^ (O^+Co) 

0 L tan r Q 

ze 2 


4 T 

s 

J mcr 


C = V 


ar,„ 

dV 


/ r d(s /fl 

da 


and 

r„ trim angle of climb, deg 

T thrust, lb 

I y pitching moment of inertia, slug-ft 2 

8 


m mass of airplane, slugs 

c wing chord, ft 

S wing area, sq ft 

u= AV/V 

V velocity, ft/sec 

AV change in velocity, ft/sec 
C L lift coefficient 

a angle of attack, radians 

C m pitching-moment coefficient 

C n drag coefficient 

p air density, slugs/cu ft 

Q 

Solving equations (Al) for - by the method of determi- 
nants gives 



x « + t T) 

■r a 

0 



z a -j- tD 

— 

8 

c m 

" l u 

c m<x \-C m[)a D 

— O 
m s 

S~' 

X u +tD 

Xa 

X 6 


z u 

Za-'rrD 

Zo—tD 


Cm u 

^m a -\-C"i Dc P 

C mM D-ZyKD* 


Expanding these determinants by the usual methods gives 
an expression for ^ consisting of the. ratio of two polyno- 
mials in D, which may be written as 
6 b\D-\- Ci 

8 u>I)* b>[)''-\ cJr | </.,!) -\ c>2 

The phase and amplitude of the response of 6 to the motion 
5=sin c of is obtained by substituting for D in the above 
0 

expression for • This substitution gives a complex number, 

say, The angle of lead of 6 ahead of 8 is tan -1 

and the amplitude of 6 is ^ A 2 +i? 2 ; that is 

, -i B 
t T ~ — tan . 

U 1 

R= y /A r +B r2 

and 

8=R sin (cd—ej) 



^ to 



APPENDIX B 


DETERMINATION OF STABLE SIDE OF STABILITY BOUNDARY OBTAINED BY FREQUENCY-RESPONSE METHOD 


It is shown in the body of this report that a critical value 
of the control gearing exists if at some frequency w the values 
of € t and Cp are equal and that the value of the critical control 
gearing is equal to 8 r j8 vi at that frequency. It will now be 
shown that the stable side of the stability boundary may be 
determined by the relative slopes of the e r and t p curves at 
the point where e r =e„. 

As shown in reference 2, the response of an airplane to any 
disturbance may be computed from the “response function.” 
The response to a sinusoidal disturbance is obtained by sub- 
stituting iu for the variable D in the response function. 
The response functions for the airplane and autopilot may 
be combined to form the stability equation of the airplane 
plus autopilot as follows: 

Let 9=Ji(D)5 bo the airplane response function and 
8 = kj 2 {D)Q be the autopilot response function where k is the 
control gearing. Combining the two gives 

8=kf 2 (D)f I (D)8 

or 

1 — kf 2 (D)f l (D) =0 

By substituting /(ZJ) =/, (D)J 2 (D) this equation may be writ- 
ten as 

\-kf(D)=Q (Bl) 

This is the “stability equation,” the roots of which for the 
variable D determine the damping and frequency of the 
motion. The frequency responses are obtained from the 
response functions by use of the following: 

Z, (ho) = Re~ Ur = \ e~ u ' 

0 r 

/ 2 (ho)=S Pi e ie » 

Then a dynamic-stability boundary exists if 1— fc/(ho) = 0 
for some value of to, say oi cr . It follows from the definition 
of /, and / 2 that this dynamic-stability boundary occurs when 

e r — e„ — 0 

and 

h=* T =k cr 

s 

at some value of co = aj er . Under these conditions, equation 
(Bl) has a pair of roots D= ±iw eT when k has the value given. 
In order to find the stable side of the boundary, it is neces- 


sary to find the sign of the real part of at the boundary. 
From equation (Bl) 

dD j(D) 
dk kj’ (D) 

This equation must be evaluated when ZZ— ico cr . At this 
value of co, f{iu> „) = }-• In general, 




Then 



If — J is negative, the real part of ^ is positive. 

As k increases above the critical value, therefore, the roots 
of the stability equation indicate instability; that is, if 



then the stable region is located where k is less than lc cr and 
vice versa. 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) ; 

Angular 

Longitudinal 

X 

X 

Rolling 

L 

Y — ->Z 

Roll 


u 

V 

Lateral 

Y 

Y 


M 

z — >x 


0 



Normal 

z 

z 


N 

x >Y 



to 

r 











Absolute coefficients of moment 

/ 1 A ri Af 

! ~qbS m ~qcS 

(rolling) (pitching) 


° n ~ZbS 

(yawing) 


Angle of set of control surface (relative to neutral 
position), 8. (Indicate surface by proper sub serin tJ 


D 

P 

P/D 

V' 

v, 

T 

Q 


Diameter 
Geometric pitch 
Pitch ratio 
Inflow velocity 
Slipstream velocity 

Thrust, absolute coefficient Or— 


4. PROPELLER SYMBOLS 
P 


P 


T 


Torque, absolute coefficient C a = 




G. 

v 

n 

$ 


prPIP 

6. NUMERICAL RELATIONS 


Power, absolute coefficient C!p= 

Speed-power coefficient = ^ p^ 2 
Efficiency 

Revolutions per second, rps 
Effective helix angle — tan^ 1 ^^^ 


1 hp=76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower =0.9863 hp 
1 mph=0.4470 mps 
\ mps=2.2369 mph 


1 lb=0.4536 kg 
1 kg =2. 2046 lb 
1 mi= 1,609.35 m= 5,280 ft 
1 m=3.2808 ft 




